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Abstract 

The large-time behavior of solutions to Burgers equation with small viscosity is de- 
scribed using invariant manifolds. In particular, a geometric explanation is provided for 
a phenomenon known as metastability, which in the present context means that solutions 
spend a very long time near the family of solutions known as diffusive N-waves before 
finally converging to a stable self-similar diffusion wave. More precisely, it is shown that 
in terms of similarity, or scaling, variables in an algebraically weighted L?' space, the self- 
similar diffusion waves correspond to a one-dimensional global center manifold of stationary 
solutions. Through each of these fixed points there exists a one-dimensional, global, at- 
tractive, invariant manifold corresponding to the diffusive N-waves. Thus, metastability 
corresponds to a fast transient in which solutions approach this "metastable" manifold of 
diffusive N-waves, followed by a slow decay along this manifold, and, finally, convergence 
to the self-similar diffusion wave. 



*einail: MargaretJ3eck@brown.edu; The majority of this work was done while MB was affiliated with the 
Department of Mathematics, University of Surrey, Guildford, GU2 7XH, UK 
''^email: cew@math.bu.edu 



1 



1 Introduction 



It is well known that viscosity plays an important role in the evolution of solutions to viscous 
conservation laws and that its presence significantly impacts the asymptotic behavior of solu- 
tions. Much work has been done to understand the relationship between solutions for zero and 
nonzero viscosity. For an overview, see, for example, |Daf05[ ILiuOOj . With regard to Burgers 
equation, one key property is the following. If = u^{x,t) denotes the solution to Burgers 
equation with viscosity fi and = u^{x,t) denotes the solution to the inviscid equation, then 
it is known that — > in an appropriate sense for any fixed t > as — *• 0. However, for 
fixed n, the large time behavior of u'^ and is quite different, and they converge to solutions 
known as diffusion waves and N-waves, respectively. Thus, the limits — > and t ^ oo are not 
interchangeable. 

Recently, a phenomenon known as metastability has been observed in Burgers equation with 
small viscosity on an unbounded domain |KT01j . Generally speaking, metastable behavior is 
when solutions exhibit long transients in which they remain close to some non-stationary state 
(or family of non-stationary states) for a very long time before converging to their asymptotic 
limit. In |KT01j . the authors observe numerically that solutions spend a very long time near a 
family of solutions known as "diffusive N-waves," before finally converging to the stable family 
of diffusion waves. This terminologjQ is due to the fact that the diffusive N-waves are close 
to inviscid N-waves. In [KTOl] this is proven in a pointwise sense. Furthermore, in terms of 
scaling, or similarity, variables, they compute an asymptotic expansion for solutions to Burgers 
equation with small viscosity. They find that the stable diffusion waves correspond to the first 
term in the expansion, whereas the diffusive N-waves correspond to taking the first two terms. 
Thus, by characterizing the metastability in terms of these diffusive N-waves, they provide a 
way of understanding the interplay between the limits /i — * and t — oo. 

In this paper, we show that the metastable behavior in the viscous Burgers equation, described 
in |KT01j . can be viewed as the approach to, and motion along, a normally attractive, invariant 
manifold in the phase space of the equation. In terms of the similarity variables, we show that one 
has the following picture. There exists a global, one-dimensional center manifold of stationary 
solutions corresponding to the self-similar diffusion waves. Through each of these fixed points 
there exists a global, one-dimensional, invariant, normally attractive manifold corresponding to 
the diffusive N-waves. For almost any initial condition, the corresponding solution of Burgers 
equation approaches one of the diffusive N-wave manifolds on a relatively fast time scale: r = 
(9(1 log Due to attractivity, the solution remains close to this manifold for all time and 
moves along it on a slower time scale, r = 0(l//i), towards the fixed point which has the same 
total mass. Note this this corresponds to an extremely long timescale t ^ 0{e^^^) in the original 
unsealed time variable. This scenario is illustrated in Figure [H 

Note that, in [LiuOOj . it was shown that the large time behavior of solutions to a general class of 
conservation laws is governed by that of solutions to Burgers equation. Roughly speaking, this is 
due to the marginality of the nonlinearity in the case of Burgers equation and the fact that any 

-'^ These diffusive N-waves are also discussed in jWhi99[ §4.5], where they are referred to simply as N-waves. 
Here, as in |KT01| . we reserve the term N-wave for solutions of the inviscid equation and diffusive N-wave for 
solutions of the viscous equation. 
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Figure 1: A schematic diagram of the invariant manifolds in the phase space of Burgers equation, 
(12.31) . and their role in the metastable behavior. The solution trajectory (blue) experiences an 
initial fast transient of t = 0{ \ log/x|) before entering a neighborhood of the manifold of diffusive 
N-waves (green). It then remains in this neighborhood for all time as it approaches, on the slower 
time scale of r = 0{l/fi), a point on the manifold of stable stationary states (red). 

higher order nonlinear terms in other conservation laws are irrelevant. Therefore, the present 
analysis for Burgers equation could potentially be used to predict and understand metastability 
in other conservation laws with small viscosity, as well. 

We remark that a similar metastable phenomenon has also been investigated in Burgers equation 
on a bounded interval [BKS951 ISW99] and numerically observed in the Navier-Stokes equations 
on a two-dimensional bounded domain |YMC03j . Furthermore, metastability has been observed 
in reaction-diffusion equations: for example, on a bounded interval [FH89t ICP89t ICP90j and 
spatially discrete lattice [G VV95j . In |CP90j . the metastable states were described in terms of 
global unstable invariant manifolds of equilibria. 

The remainder of the paper is organized as follows. In ^ we state the equations and function 
spaces within which we will work, as well as some preliminary facts about the existence of 
invariant manifolds in the phase space of Burgers equation. We also precisely formulate our 
results. Theorems [T] and [21 In sections §31 SI and ^ we prove Lemma 12.11 Theorem [H and 
Theorem [21 respectively. Concluding remarks are contained in ^ Finally, the appendix contains 
a calculation that is referred to in ^ 



2 Set-up and statement of results 

We now explain the set-up for the analysis and some preliminary results on invariant manifolds. 
Our main results are precisely stated in §2.31 
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2.1 Equations and scaling variables 



The scalar, viscous Burgers equation is the initial value problem 

dtu = fid^u — uUx 
u\t=o = h, (2.1) 

and we assume the viscosity coefficient fi is small: < ;U ^ 1. For reasons described below, it 
is convenient to work in so-called similarity or scaling variables, defined as 

« = r = log(« + 1). (2.2) 

In terms of these variables, equation (12.11) becomes 

drW = C^w — ww^ 
w\r=o = h, (2.3) 

where C^w = fid'^w + \d^{^w). 

We will study the evolution of (12.31) in the algebraically weighted Hilbert space 
L\m) = 1^ G L\R) : = J {1 + eTHOl'd^ < oo} . 

It was shown in |GW02j that, in the spaces L'^{m) with m > 1/2, the operator generates a 
strongly continuous semigroup and its spectrum is given by 

a{C,) = G n} U |a G C : Re(A) <\-j]- (2-4) 

This is exactly the reason why the similarity variables are so useful. Equation (12.41) shows that 
the operator £^ has a gap, at least for m > 1/2, between the continuous part of the spectrum 
and the zero eigenvalue. As m is increased, more isolated eigenvalues are revealed, allowing 
one to construct the associated invariant manifolds (see below for more details). In contrast, 
the linear operator in equation (12. ip . in terms of the original variable x, has spectrum given by 
(—00,0], which prevents the use of standard methods for constructing invariant manifolds. 

For future reference, we remark that the eigenfunctions associated to the isolated eigenvalues 
A = —n/2 are given by 

cp,(^) = -^e~i, ^„(0 = (5F^o)(0- (2.5) 

See |GW02] for more details. Smoothness and well-posedness of equations (12. ip and (12. 3p can 
be dealt with using standard methods, for example information about the linear semigroups and 
nonlinear estimates using variation of constants. 
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2.2 Invariant manifolds 



We now present the construction, in the phase space of (12.31) . of the exphcit, global, one- 
dimensional center manifold that consists of self-similar stationary solutions. We remark that 
this is similar to the global manifold of stationary vortex solutions of the two-dimensional Navier- 
Stokes equations, analyzed in [GW05j . First, note that stationary solutions satisfy 

ddf^w^ + ^^w - ^w"^) = 0. 
They can be found explicitly by integrating the above equation and rewriting it as 

Integrating both sides of this equation from — oo to ^ leads to the following self-similar stationary 
solution, for each G M: 



1_ ajj^e-^V(4/.)(i^ 

Note that equation (12.31) preserves mass, and we can therefore characterize these solutions by 
relating the parameter to the total mass M of the solution. We have 

M= / w{i)di = a^ / —- TJ-r—di =-2/ilog(l-«oW-), 

where we have made the change of variables 6* = 1 — |^ /f^ e'^'^/^^'^-'cir/. Therefore, we define 



These solutions are often referred to as diffusion waves [LiuOOj . 

For m > 1/2, the operator has a spectral gap in L'^{m). By applying, for example, the results 
of |CHT97j . we can conclude that there exists a local, one-dimensional center manifold near the 
origin. In addition, because each member of the family of diffusion waves is a fixed point for 
(12.31) . they must be contained in this center manifold. Thus, this manifold is in fact global, as 
indicated by figure [TJ 

Remark 2.1. Another way to identify this family of asymptotic states is by means of the Cole- 
Hopf transformation, which works for the rescaled form of Burgers equation as well as for the 
original form (12.11) . If w is a solution of (12. 3p . define 

W{^,t) = w{^,T)e^^J--'"^y'^^'y = -2ftd^exp (-^ w{y,T)dy^ . (2.7) 

A straightforward computation shows that W satisfies the linear equation 

drW = L^W. (2.8) 
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Conversely, let W be a solution of (12. 8p for which 1 — 2^ W{y, T)dy > for all ^ and 
T > 0. Then the inverse of the above Cole-Hopf transformation is 



w{i, t) = -2/i95 log (1 - ^ / wiy^ ^)dy) ■ 



(2.9) 



The family of scalar multiples of the zero eig en function, /3oV^o(0> where ip^ is given in l{2.5\) . is 
an invariant manifold (in fact, an invariant subspace) of fixed points for (\2.8\\ . Thus, the image 
of this family under (12. 9p must be an invariant manifold of fixed points for (12.31) . Computing 
this image leads exactly to the family \2. 0\] . where l3o = V^vr/Iao. 

Remark 2.2. One can prove that this self-similar family of diffusion waves is globally stable 
using the entropy functional 



H[w]{t)= I w{^,T)e-^^-^'"^y'^^''Ho, 



This is just the standard Entropy functional for the linear equation ^2. <SI) with potential (4/i), 
in combination with the Cole-Hopf transformation. For further details regarding these facts, see 
IDiFOSf . 

We next construct the manifold of diffusive N- waves. Recall that, by (12.41) . if m > 3/2 then 
both the eigenvalue at and the eigenvalue at —1/2 are isolated. The latter will lead to a one 
dimensional stable manifold at each stationary solution. 



To see this, define w = Am + v and obtain 



Affv - vv^ 



(2.10) 



where is just the linearization of (12.31) about the diffusion wave with mass M. One can 
see explicitly, using the Cole-Hopf transformation, that the operators and A'^^ are conjugate 
with conjugacy operator given explicitly by 



U = d^ 



Af^'U = UC, 



(2.11) 



Thus, one can check that the spectra of the operators Aff and £^ are equivalent in L'^{m). 
Furthermore, we can see explicitly that the eigenfunctions of A'^ are given explicitly by 



$n(0 = 



/-oo Vn{.y)dy 



where <y9„ is an eigenfunction of Notice that, up to a scalar multiple, $1 = O^Am- If we 
chose m > 3/2, by (12.41) we can then construct a local, two-dimensional center-stable manifold 
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near each diffusion wave. We wish to show that, if the mass is chosen appropriately, then this 
manifold is actually one-dimensional. Furthermore, we must show that this manifold is a global 
manifold. 

To do this, we appeal to the Cole-Hopf transformation. Using (12.71) . we define 

and find that V solves the linear equation 

d,V = AffV. 

Thus, the two-dimensional center-stable subspace is given by span{$07 The adjoint eigen- 
function associated with $o is just a constant. Therefore, if we restrict to initial conditions that 
satisfy 

V(e,OK = 0, (2.12) 



then the subspace will be one dimensional and given by solutions of the form 

y(e,r) = ai$i(Oe-i 
One can check that condition (I2.12p is equivalent to 

Jr 

Since w = Am + v, we can insure this condition is satisfied by choosing the diffusion wave that 
satisfies 



M= f AM{Od^= [ w{^,0)d^. 
Jr Jr 



Thus, near each diffusion wave of mass M, there exists a local invariant foliation of solutions 
with the same mass M that decay to the diffusion wave at rate e~2'^. 

To extend this to a global foliation, we simply apply the inverse Cole-Hopf transformation (12.90 . 
as in Remark 12.11 to the invariant subspace 

{V{^,t) = ai<l>i(0} = {V{^,t) = a^^^AM}. 
This leads to the global stable invariant foliation consisting of solutions to (12.101) of the form 



Using the relationship between v and w, this foliation leads to a family of solutions of (12.31) of 
the form 

WNi^, r) = A.siO + M^, r) = A^iO + , "'!73'^^Lv (2-13) 

Below it will be convenient to use a slightly different formulation of this family, which we now 
present. 
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The subspace spanjy^o, V^i}, corresponding to the first two eigenfunctions in fl2.5p . is invariant 
for equation fl2.8p . Therefore, as in Remark 12. 11 by using the inverse Cole-Hopf transformation 
(12.91) we immediately obtain the explicit, two parameter family 

u.„(f,r) = ftv'o(0+ft^-^'^i(0 (2,14) 

where 



Po = V^ao = 2fi{l-e--^). (2.15) 

Based on the above analysis, (I2.14p and (I2.13P are equivalent. Note that, although the method 
used to produce (I2.14p is much more direct than that of (I2.13p . we needed to use the operator 
Am and its spectral properties to justify the claim that this family does in fact correspond to 
an invariant stable foliation of the manifold of diffusion waves. 

We now explain why solutions of the form (12.140 are referred to as the family of diffusive N- waves. 
As mentioned in §1, this terminology was justified in |KT01j by showing that each solution wn 
is close to an inviscid N-wave pointwise in space. Since we are working in L'^{m), we need to 
prove a similar result in that space. 

Recall some facts about the N-waves, which can be found, for example, in |LiuOO] . Define 

/y r'xi 
u{x)dx, and g = 2sup / u{x)dx, (2-16) 
-oo y Jy 

which are invariant for solutions of equation ( 12. ID when /i = 0. (Note that our definitions of p 
and q differ from those in [KTOl] by a factor of 2.) The mass satisfies M = {q — p)/2. We will 
refer to q as the "positive mass" of the solution and p as the "negative mass" of the solution. 
The associated N-wave is given by 



Xp,q{x,t) 



if - ^p{t + 1) <x< ^qit + 1) 



t+i 

otherwise 



which is a weak solution of (12. ip only when /i = 0. When < /i <C 1 it is only an approximate 
solution because the necessary jump condition associated with weak solutions is not satisfied. 
One can check that its positive and negative mass are given by q and p. In terms of the similarity 
variables (12.20 . this gives a two-parameter family of stationary solutions 



NUO = \l 'f-v^<«<v^ (2.17) 

I otherwise 

of equation (12. 3p when /i = 0. 

We now relate the quantities (5q and (3i in (I2.14p to the quantities p and q. These calculations 
follow closely those of [KTOlt §5]. Using equation (I2.15P and the fact that M = {q — p)/2, we 
see that 

/5o = 2Ml-e" — ) = <^ ^ ' ^ (2.18) 

-2/ie 4" + 0{fi) if q <p, 
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where exp = 0{e '"^^) for some C > 0. Using the calculation in the appendix, one can relate 
the quantity Pi in fl2.14p . for any fixed r, to the quantities p and q via 

/?,e-^ = -V/2 0FeP/(4/^) - 1^ + for < g < p, (2.19) 

and a similar result holds for g > p > 0. Two key consequences of this, which will be used 
below, are 



The quantities /?o and /3i are related via 

^ = exp, (2.20) 

Pi 

The values of p and q for the diffusive N- waves change on a timescale of r = 0(^). (Recall 
they are only invariant for /i = 0.) 



This second property, which can be seen by differentiating f l2.19p with respect to r, will lead to 
the slow drift along the manifold of diffusive N-waves (see below for more details). 

The following lemma, which will be proved in ^ states precisely that there exists an N-wave 
that is close in L'^{m) to each member of the family w^, at least if the viscosity is sufficiently 
small, thus justifying the terminology "diffusive N-wave." 

Lemma 2.1. Given any positive constants 6, p and q, let wn{^iT) be a member of the family 
Ili2.14\ ) of diffusive N-waves such that, at time r = Tq, the positive mass ofwi\[{-, Tq) is q and the 



negative mass is p. There exists a > sufficiently small such that, ifO<n< hq, then 

\\wN{-,ro) - Np,g{-)\\L2(m) <S 



2.3 Statement of main results 



We have seen above that the phase space of fl2.3p does possess the global invariant manifold 
structure that is indicated in figure [TJ To complete the analysis, we must prove our two main 
results, which provide the fast time scale on which solutions approach the family of diffusive N- 
waves and the slow time scale on which solutions decay, along the metastable family of diffusive 
N-waves, to the stationary diffusion wave. 

Theorem 1. (The Initial Transient) Fix m > 3/2. Let w{^,t) denote the solution to the 
initial value problem Ii2.3\) whose initial data has mass M, and let iVp^q(^) be the inviscid N-wave 
with values p and q determined by the initial data w{^,0) = h{^) € L?'{m). Given any 5 > 0, 
there exists aT > 0, which is 0{\ log/i|), and /i sufficiently small so that 

\\w{T)-NpJ\l^„,)<S. (2.21) 
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This theorem states that, although the quantities p and q are determined by the initial data 
w{^, 0), w is close to the associated N-wave, iVp g, at a time t = T = 0{\ log/i|). The reason for 
this is that p = p{t) and q = q{T) change on a time scale of 0{l/fi), which can be seen using 
equation (I2.19P and is slower than the initial evolution of w. The rate of change of p and q also 
determines the rate of motion of solutions along the manifold of diffusive N-waves, as illustrated 
in Figure [H Note that this theorem states that the solution will be close to an inviscid N-wave 
after a time T = 0{ \ log/i|). By combining this with Lemma [2711 we see that the solution is also 
close to a diffusive N-wave. 

We remark that the timescale (9(1 log was rather unexpected. We actually expected to 
approach a diffusive A^-wave on a time scale 0{1), although we have not yet been able to obtain 
this stronger result. However, these time scales correspond well with the numerical observations 
of [KTOll Figure 1], where one can see that, for fi = 0.01, the solution looks like a diffusive 
N-wave at time 2 and a diffusion wave at time 100. 

Remark 2.3. To some extent, this fast approach to the manifold of N-waves can be thought of 
in terms of the Cole- Hopf transformation 1^2. 7\ ), which depends on fi. For small fi, this nonlinear 
coordinate change can reduce the variation in the solution for \^\ large. This is illustrated in 
figure 5. 1 of IKNO0^ . If fi is small enough, the Cole Hopf transformation can make the initial 
data look like an N-wave even before any evolution has taken place. 

Theorem 2. (Local Attractivity) There exists a cq sufficiently small such that, for any 
solution w{-,t) of the viscous Burger's equation 1^2. 3\} for which the initial conditions satisfy 

W\r=0 = W% + 0°, 

where w% is a diffusive N-wave and \\(f)'^\\L^(m) < co, there exists a constant such that 

w(-,r) = WAr(-,r) + 0(-,r), 
with wn the corresponding diffusive N-wave solution and 

||0(-,7-)||i2(„) < C^e~^. 

By combining these results, we obtain a geometric description of metastability. Theorem [1] and 
Lemma [2?T] tell us that, for any solution, there exists aT = 0{\ log at which point the solution 
is near a diffusive N-wave. By using Theorem [2] with this solution at time T as the "initial 
data," we see that the solution must remain near the family of diffusive N-waves for all time. As 
remarked above, the time scale of 0{l/n) on which the solution decays to the stationary diffusion 
wave is then determined by the rates of change of p{t) and g(r) within the family of diffusive 
N-waves. In other words, near the manifold of diffusive N-waves, w{^,t) = wiy{^,T) + 0(^,r), 
where r) ~ e~'^ and WAr(^,r) is approaching a self-similar diffusion wave on a timescale 
determined by the rates of change of p(r) and qi^r), which are 0{l/fi). 

We remark that it is not the spectrum of that determines, with respect to fi, the rate of 
metastable motion. Instead, this is given by the sizes of the coefficients /?o and /5i in the spectral 
expansion and their relationship to the quantities p and q. 
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3 Proof of Lemma 



2.1 



We now prove Lemma I2.1[ 

In |KT01] . Kim and Tzavaras prove that the inviscid A^-wave is the point-wise hmit, as ^ 0, 
of the diffusive iV-wave. Here we extend their argument to show that one also has convergence 
in the L^{m) norm. For simphcity we will check explicitly the case in which < q < p - the case 
in which q is larger than p follows in an analogous fashion. However, we note that we do require 
that both p and q be nonzero which is why we stated in the introduction that our results hold 
only for "almost all" initial conditions. See Remark 13.11 below. 

Using equation (12.131) we can write the diffusive A^-wave with positive and negative mass given 
by q and p at time tq as 

, PoMO + /5ie-^o/Vi(0 PoMO + PmiO 



1 - I /-oo My)dy - |e--/Vo(0 1 - I IL My)dy - |v^o(0 

where for notational simplicity we have defined Pi = Pie~'^°^'^. If we now recall that v'llO — 
— ^(/?o(05 ^6 can rewrite the expression for the as 

yJN{tT„) = - „. — — . (3^22) 



We need to prove that 



We'll give the details for 



(l + eT(^^(e,ro)-iV,,,(Or 



The integral over the negative half axis is entirely analogous. 

Break the integral over the positive axis into three pieces - the integral from [0, ^/q — e], the 
integral from [^/q — e, ^/g + e] and the integral from [y/q + e, 00). Here e is a small constant that 
will be fixed in the discussion below. We refer to the integrals over each of these subintervals as 
/, //, and III respectively and bound each of them in turn. 

The simplest one to bound is the integral II. Note that, using equations fl2.15p and (12.191) . the 
denominator in (I3.22p can be bounded from below by 1/2 and, thus, the integrand is can be 
bounded by C(l + {^/q + e)^)'"^. Therefore, if e < ^^q, we have the elementary bound 

// < Ceg(l + 4g)"\ 
Next consider term ///. For ^ > y/q, Np^q{^) = so 

POO 

111= / {i+er{w,,,{^,To)Yd^ (3.23) 
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To estimate this term we begin by considering the denominator in (13.221) . Using (12.181) and 
(12191) . we have 

1 _^ /"^ 



My)dy = e-^'"'"' + ^ / ipo{y)dy. (3.24) 
Thus, the full denominator in (13.221) has the form 

2/^ J 



- 1 + ^v^^'^' ei(P-'')e«'/(^'^) + exD 

^ 1 v/^/i3/2eP/{4M) + O(^) ^ ^""P 



1 + C(^-ie-p/(4A')) 



exp, 



where exp denotes terms that are exponentially small in /i (i.e. contain terms of the form 
g-p/(4M) QY g-'j/C^M))^ uniformly in ^. Note that, in the above, the term (fQ{y)dy/ipo{^) was 
bounded uniformly in /x using the estimate 



^2 I ^2 

e~~ds < —e~~, for x > 0, 

X 



which can be found in |KS91l Problem 9.22]. But with this estimate on the denominator of t^jv, 
we can bound the integral /// by 



oo 



iii<c I (1 + en^ - ^)'(i + e^(«'-''))-2rfe, 

where the constant C can be chosen independent of fi for /i < /io if yUo is sufficiently small. This 
integral can now be bounded by elementary estimates, and we find 

/// < Ce-'^/^", 

where the constant C depends on q but can be chosen independent of /i. 
Finally, we bound the integral /. Note that for < ,^ < y^g, Np^q{^) = ^, so 



Wp,qi^, To) - iVp,g(0 



.2^+^e4k«'-«)+^exp 



(1 + e(«'-'?)/4/^ + exp) 

However, using our expressions for (3o,i in terms of p and q and the fact that in term I C,^ — q < 
—e^/q we see that all of these terms are exponentially small. Since the length of the interval of 
integration is bounded by ^ we have the bound 

+ exp. 



Combining the estimates on the terms /, //, and /// we see that if we first choose e << 6 and 
then fi << e, the estimate in the lemma follows. This completes the proof of Lemma [2.11 

Remark 3.1. The calculation in the appendix shows that Pi = if and only if p = or q = 0. 

Therefore, in that case, is really just a diffusion wave, and so there is no metastable period 
in which it looks like an inviscid N-wave. 
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4 Proof of Theorem [T] 



In this section we show that for arbitrary initial data in L^(m), m > 3/2, the solution approaches 
an inviscid N-wave in a time of 0{ \ log/i|), thus proving Theorem [H 

Remark 4.1. Here we will use a different form of the Cole-Hopf transformation than that given 
in (2.1). In particular, we will use 



f/(x,t) =e-^^— (4.25) 



Equation (2.1) is essentially the derivative of lj^4-25\ ), and both transform the nonlinear Burgers 



equation into the linear heat equation. Each is useful, for us, in different ways. Equation ^2.7\ ) 
preserves the localization of functions, for example, whereas ^.25^ leads to a slightly simpler 
inverse, which will be easier to work with in the current section. 



Using the Cole-Hopf transformation (14.251) and the formula for the solution of the heat equation, 
we find that the solution of (12.11) can be written as 

U(x,t) = — r , 

where H{y) = h{z)dz. If we change to the rescaled variables (12. 2p . this gives the solution 
to (12. 3p in the form 

w ^ = — ^TTi -■ 4.26 

^ Je~U2(^-^^'+''^^^^'^^)dr] 



We will prove that, if we fix 5 > 0, then there exists a fi sufficiently small and a T sufficiently 
large {0{ \ log/i|) as /i — ^ 0) such that 

||W(-,T) - Np^qi-)\\L^m) < 

We estimate the norm by breaking the corresponding integral into subintegrals using (— oo, —^/p— 

e), {-Vp - ^' -Vp + (~Vp + ^' ~^)' (~^' '^)' - e), - e, + e) and (^ + e, cx)). 

Note that the integrals over the "short" intervals can all be bounded by Ce, so we'll ignore 
them. We'll estimate the integrals over (e, ^/q — e) and (y^+ ^) oo) - the remaining two are very 
similar. 

First, consider the region where ^ > ^/q + e. In this region, N{C,) = 0, so we only need to show 
that, given any 6 > 0, there exists a ^ sufficiently small and T > 0, of 0{\ log/i|), such that 

/■oo 

/ {i+erH^,r)\'d^<s. 

Consider the formula for w, given in equation (I4.26p . To bound this, we must bound the 
denominator from below and the numerator from above. We will first focus on the denominator. 
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We will write the denominator as 



CXD 

/*oo 

= Ji + /2 + /a, 

for some i? > that will be determined later. Consider the first integral, Ji. In this region 

oo (i + r)2 

1 (4.27) 



oo 



where the constant C(-R, ||/i||m) — * as i? ^ oo or ||/i||m — ^ 0. In addition, note that the error 
function satisfies the bounds 

e 2 < / e 2 rfs < _e 2 



for z > [KS91t pg 112, Problem 9.22]. Therefore, we have that 

_ C{R,m\m) f~^"^ „_L(£_„)2 

/i > e / e '''' rf?7 

J — oo 

-c(R,\\h\\m) 2u(f + i?e"^/^) (e+-Re~i)2 
2/i+(e + i?e-2)2 

^ U -C(R,\\h\U) (S + fle~i)2 

> C—e 2m e 4m 



(4.28) 



— ( 



^ /i -C(fl,||fe||m) fl2^-r 

> C — e 2m e 2m e 2m. 



where we have used the fact that (a + 6)^ < 2a^ + 26^. In order to bound I2, we will use that, 
for T] > Re2, similar to f l4.27p 



Therefore, 



e^rj POO 

h{y)dy\ = \M-l h{y)dy\<M + CiR, 

■ 00 J e'^rj 



poo 



_,2 , (4.29) 
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Note that in making the above estimate, we have chosen r large enough so that \Re a] < e/2, 
and so ^ - Re'^^^ > 0. Thus, the error function is bounded from below by -^/vr/^- 



Consider h. We can bound 



Therefore, 



H{eir]) = M-l h{y)dy<M+-. 



/"^""^ , (4.30) 



Taking the largest of equations (gSSD, (1122]), and (jl^ . we obtain 

y"e-27i(5(5-'?)'+^^(-^/'^)) > Cv/A^e-57i(^+^(«'ll''ll™)). (4.31) 

Now, we will bound the numerator in equation (14.260 from above. We will split the integral up 
into the same three regions as above, and denote the resulting terms by Ji, J2, and J3. First, 
we have 

I Jil = I / (e - r/)e-^(«-'')'e-i^(^^'')rfr7| 



< |e 2m / (^ - ry)e ^^^^ 



00 

2 



(4.32) 



= e 2m 2/ie 

C(R,\\h\\m) 

< 2/ie 2m e 4m. 



Next, consider J2. We have 



I J2I = I / (^ - r/)e-^(«-^)'e-^^(^^'')c/r/| 



Re" 7 



If we now integrate by parts inside the integral and use the fact that if(e2r/) > M — {q/2), we 
obtain 

/•oo ^ 

I J2I < C/ie-iei[^-(«-^^"^'')'] + e-i| / e-i(«-'')'/i(e-/2r/)e+i(*^-^('=^''))rfr/|. (4.33) 

We now turn to J3. Again we use the fact that H{e^r]) > M — {q/2). Then 

I J3I < 2Ci?e-5e-i(e + i?e-5)e^[''-(«-^'^"^'')'l. (4.34) 
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Combining the estimates for the Jj's, equations (14.321) - fl4.34p . and the estimate for the denom- 
inator (14.311) . we have 

/•OO noo 2 

/ {l + er\w{^,r)\'d^<C / (l + e')>e^ef^(^'ll"ll™)e-fcrfe 



OO 



/•OO 

+ C / (1 + e2)>e^^(^'ll'^ll™)e^[^-(«-^^"^'')']rfe 

/•OO -| /'OO ^ 

= / + // + /// + JV. 

We now estimate term II. Terms I and III are similar. Define z = ^ ~ ^ ^ (0, oo). Recalling 
that r has been chosen sufficiently large so that Re~^ < e/2, we have 



Therefore, we have 
C 



io 2 



This can be made as small as we like (for any q) by choosing R large enough so that C {R, \ \h\ \m) < 
e^/16 and fi sufficiently small. 



Term IV can be bound by 



|jy|<C-et^(^'ll'^ll-)||/(.)*^?(z)||i„ 



where f{z) = (1 + \z\'^)h{e'^'''^z) and g{z) = (1 + \z\'^)e~^^ , and we have used the fact that 
(1 + 1^1™) < C(l + 1^ — ri\"^){l + |?7|™). Estimating the convolution by the norm of g and the 
norm of /, we arrive at 

|/\/| <Ce-^/2gfc(ii,||h|U)||/^||2^_ 

In order to make this term small, we much chose R so that C{R, \ \h\\m) ~ /i as /i — > 0. One 
can check that C{R, \ \h\\m) < C\\h\\m/ R"^"^'^- Since we have required that Re~'^^'^ < e/2, this 
means we must chose r large enough so that 

r>^|log(/.)|. 

Term IV will then be small because e""^/^ is. 
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Next consider the part of the integral contributing to ||tf(-,r) — A^(-)||m for e < ,^ < ^ — e. 
We assume, as above, that Re-^/^ < e/2. From flCTj) and the fact that N{i) = { for ^ in this 
range, we have 

As above we will split the denominator up into three pieces, Ji - J3, to bound it from below, and 
split the numerator up into three pieces, Ji - J3, to bound it from above. Many of the estimates 
are similar to those above, and so we omit some of the details. 

For the denominator, we have 

■00 

Also, 

where, as above, we have used the fact that ^ — Re~'^^'^ > 0. Finally, we have 



For the numerator, we have 



I Ji I < e 2f / -r]e ^« "i' drj 



■00 

00 2 



g(«,l|fa||m) _^ 

< e 2m / — t,)e *i^az 

ic+iJe-^/2 

C(R,\\h\\m) (e+-Re-^/2)2 (g+fl<;-^/2)2 

< C/ie 2m e 4m -I- C^^e ^m 



where we have used the fact that ^ > 0. Next, 



_M C(R,\\h\\,n) I /" 73m I 2 

J2I < e 2Me 2m / - Y/4/i2;)e ^ rf^ 

J —00 



_M C{R,\\h\\m) (g-fle-^/2)2 

< C/i.^^ 2Me 2m e 



Lastly, 

I J3I < Ce"^e "' '2m""" (i?e-^/^)^e"37^ 
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Therefore, we have 



[i+erM^,r)-N{o\'d^ 



Jo 

where we have used the change of variables ^ = z + e. This can be made small by first choosing 
R large enough so that 2C{R, \ \h\\m) < and then taking fj, small. 



5 Proof of Theorem [2 



In the previous section we saw that in a time r = 0{\ log/i|) we end up in an arbitrarily small 
(but 0{1) with respect to fi) neighborhood of the manifold of diffusive A^-waves. In the present 
section we show that this manifold is locally attractive by proving Theorem [21 

Remark 5.1. An additional consequence of the proof of this theorem is that the manifold of 
diffusive N-waves is attracting in a Lyapunov sense because the rate of approach to it, 0{e~^), 
is faster than the decay along it, 0{e~'^l'^\ Note that this is not immediate from just spectral 
considerations since this manifold does not consist of fixed points. Therefore, the eigendirections 
at each point on the manifold can change as the solution moves along it. 

Remark 5.2. In \KNO^ %5], the authors make a numerical study of the metastatic asymptotics 
of Burgers equation. Their numerics indicate that, while the rate of convergence toward the 
diffusive N-wave (e~'^ in our formulation) seems to be optimal, the constant in front of the 
convergence rate (Cfp in our formulation) can be very large for some initial conditions. In fact, 
our proof indicates that this constant can be as large as max{l, e^"'^/^'^}. 

To prove the theorem note that by the Cole-Hopf transformation we know that if w{^,t) = 
u^n{^,t) + 4>{C,,t) solves the rescaled Burger's equation, where wn is given in equation (12.141) . 
then 

W{^,t) = (i/;^(e,r) + 0(^,r))e^^^--("^(^'")+^(^'"))''^ 
is a solution of the (rescaled) heat equation: 

drW = C^W. (5.35) 
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We now write W = Vn + '^, where Vn = WNexp{-j^ J^^WN{y,T)dy) = (3o(po{0+Pie ^''Vi(0- 
That is, Vn is the heat equation representation of the diffusive iV-wave, which we know is a 
hnear combination of the Gaussian, ip^, and ipi. 

With the aid of the Cole-Hopf transformation we can show that decreases with the rate 
claimed in the Proposition. To see this, first note that if we choose the coefficients and f5i 
appropriately we can insure that 



vi/(e,o) = J e^(e,o)de = o. 



This follows from the fact that the adjoint eigenfunctions corresponding to the eigenfunctions 
ipo and (fi, respectively, are just 1 and — ^. This in turn means that there exists a constant 
such that 

||^(-,r)|U2(^) <C^e-^ (5.36) 
at least if m > 5/2. Integrating the Cole-Hopf transformation we find 

^ (Viv(y,r) + ^(y,r))rfy=-2/i|e^57i/-».«"'^fe'-)+*(s''-))'^2'_l) (5.37) 



and, in the case = 0, jf^ = — 2/i |e I-oo'^n ^ "'"}■ later use we note the following 
easy consequence of fl5.37p : 

Lemma 5.1. There exists a constant 6^ > such that for all t > we have 

— {VN{y, t) + ^{y, r)) dy = e"^ f-^My,r)+Hy,r))dy > _ 
-oo 2/i 

Proof. For any finite r the estimate follows immediately because of the exponential. The only 
thing we have to check is that the right hand side does not tend to zero as r ^ oo. However 
this follows from the fact that we know (from a Lyapunov function argument, for example) that 
wn{^, t) + 4>{C,, t) — > ^m(0 as r — *• cxd, where Am is one of the self-similar solutions constructed 
in Section [1] and hence 



□ 



Next note that by rearranging fl5.37p we find 



(f){y,T)dy = -2fi\og< ^ (5.38) 

I ^-2tlJ-ooVN{y,r)dy J 

Differentiating, we obtain the corresponding formula for 0, namely, 

, 1 *(e,r)/i^V;v(l/,r)rfy-V^(e,r)/«^vl/(^,r)rfy-2/iVl/(e,r) 
(p[q,T) = y ^-p r- (5.39) 

2^ (l - 4 lUVNiy, r) + M/(y, r))dy) (l - ^ /f ^ ^^(y, r)dy) 
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But now, by Lemma [STTl the denominator of the expression for can be bounded from below by 
while in the numerator Vn and Vn are bounded in time while ^ ^ and ^ are each 
bounded by C^^^e'"^ which leads to the bound asserted in Theorem [2l 



6 Concluding remarks 



In the above analysis, the global, one-dimensional manifold of fixed points governing the long- 
time asymptotics is constructed directly, in a manner similar to that of [GW05] for the Navier- 
Stokes equations in two- dimensions. However, we then utilized the Cole-Hopf transformation to 
extend the stable foliation of this manifold to a global foliation, which consists of the diffusive 
N-waves that govern the metastable behavior. Ultimately one would like to obtain a similar 
geometric description of, for example, the numerically observed metastability in [YMC03] for 
the two-dimensional Navier-Stokes equations. In order to do this one would need an alternative 
way to construct a global foliation of the manifold of fixed points. This would be an interesting 
direction for future work. 



7 Acknowledgements 



The authors wish to thank Govind Menon for bringing to their attention the paper [KTOl] , which 
lead to their interest in this problem. The second author also wishes to thank Andy Bernoff for 
interesting discussions and suggestions about this work. Margaret Beck was supported in part 
by NSF grant number DMS-0602891. The work of C. Eugene Wayne was supported in part by 
NSF grant number DMS-0405724. Any findings, conclusions, opinions, or recommendations are 
those of the authors, and do not necessarily reflect the views of the NSF. 



8 Appendix 



We now give the calculation that leads to (I2.19p . Using the definition of p in (I2.16p . we find 
that 



W7v(^,r)rf^ = 4/isup / d^log 1-—/ (po{y)dy - —e-^ipoi^) 

'OO y J -oo L ^/^ J -oo '^f^ 



= Afi sup log 
y 

A direct calculation shows that the supremum is achieved at 

2/i/3o 



y 



Pie- 



20 



Substituting in this value and rearranging terms, we find 



(po{z)dz 



Since ipQ{z)dz E (0, 1), the value of /3o in (12.151) implies that the right hand side satisfies 



if M > 0, and 



— 1 < e-if* 



p _ M_ jp_ 

g4n _ g 2m < g4M 



1 - / fo{z)dz ) < e^M - e 



if M < 0. Because (fo{y) ^ and p > 2|M|, we see that, if fi is sufficiently small, then (3i < 0. 
Furthermore, /3i = if p = — M/2, ie M < and q = 0. Using the fact that (p{y) < l/^/47^JI, 
we see that 

g4M — 1 — C 



2/i 



where c G (0, 1). This leads to the estimate fl2.19p . at least when both g 7^ and p 7^ 0. 
We remark that /3i < and the fact that 

/3o 



1 



2/i 



(po{z)dz = 1 — (1 — e 



2m 



(po{z)dz G 



(e~^,l) ifM>0 
;i,e"i) ifM<0 



implies that the denominator in the definition of wat (12.141) is never zero. 
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